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1.  INTRODUCTION 


The  aerodynamic  forces  and  moments  on  a  missile  are  produced  by  its  relative  motion  with 
respect  to  the  air  and  depend  on  the  orientation  of  the  missile  with  respect  to  the  airflow.  In 
uniform  airflow  these  forces  and  moments  are  unchanged  after  a  rotation  around  the  free-stream 
velocity  vector.  Therefore,  two  orientation  angles,  with  respect  to  the  relative  wind,  [Vg"^],  are 
needed  to  specify  the  aerodynamic  forces  and  moments.^  There  are  two  methods  of  defining  the 
two  orientation  angles.  In  aircraft  and  missile  aerodynamics,  the  Cartesian  aero-angles  (CA)  are 
most  commonly  used.  An  alternate  used  in  missile  aerodynamics  is  the  polar  aeroballistic  angles 
(PA).  Converting  CA  orientation  angles  to  PA  and  vice-versa  is  simple  but  since  the  equations  of 
motion  and  thus  force  and  moment  coefficients  associated  with  each  are  different,  the  complete 
conversion  process  is  more  complex. 

The  most  frequently  used  method  of  predicting  the  aerodynamic  coefficients  for  missiles  is  the 
Missile  DATCOM  (MD)  semi-empirical  aerodynamic  prediction  code.  In  MD  the  body 
orientation  can  be  input  in  either  of  the  two  methods  discussed  above,  but  the  aerodynamic 
coefficients  output  are  those  associated  with  CA. 

The  CA  system  works  well  with  the  equations  of  motion  for  3  degree  of  freedom  (DOF)  flyout 
simulations.  For  a  six  DOF  simulation  of  a  tetragonally  symmetrical  missile,  the  equations  of 
motion  are  simplified  if  the  PA  system  is  used.  A  program  which  converts  MD  output  from  CA 
coefficients  to  PA  coefficients  would  streamline  the  work  required  to  run  such  a  simulation  when 
wind  tunnel  aerodynamic  data  are  not  available. 


2.  REFERENCE  SYSTEMS  &  EQUATIONS 

In  converting  from  one  reference  system  to  another,  the  basis  for  each  system  must  be  defined. 
Sections  2.1  through  2.3  review  the  reference  systems,  force  and  moment  coefficients  as  well  as 
the  equations  of  motion  for  the  PA  system. 

2.1  BODY-AXIS  SYSTEM 

In  the  body-axis  system  shown  in  Figure  1,  the  axes  are  fixed  to  the  body  with  the  origin  at  the 
center  of  mass  of  the  vehicle  with  the  body  x-axis  (1®)  pointing  forward  and  the  body  z-axis  (3®) 
pointing  down.  These  axes  usually  coincide  with  the  principal  axes  of  inertia. 

2.2  CA  SYSTEM 

The  CA  orientation  angles  are  also  shown  in  Figure  1.  The  orientation  of  the  vehicle  relative  to 
the  resultant  tvind,  [Vg^],  is  defined  by  the  angle  of  attack,  a,  and  the  sideslip  angle,  (3.  The  first 
rotation  defines  the  stability  axis,  1^,  where  a  is  the  angle  between  the  body-fixed  x-axis  and  the 
stability  x-axis.  Alpha  is  positive  if  the  rotation  about  the  body-fixed  y-axis  (2^)  is  negative,  thus  a 
positive  a  is  shown.  The  second  rotation  leads  to  the  wind  axes,  where  the  sideslip  angle,  p,  is  the 
angle  between  the  stability  x-axis  (1^)  and  the  wind  axis.  Beta  is  positive  if  the  rotation  about  the 
stability  z-axis  (3^)  is  positive,  thus  a  positive  P  is  shown. 

The  force  and  moment  coefficients  for  the  CA  system  along  with  their  MD  variable  names  are 
listed  in  the  Table  1.  All  but  is  calculated  by  MD. 
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Table  1.  CA  Coefficients  and  MD  Variable  Names. 


Force  Coefficients: 

Definition: 

CA 

Axial 

CN 

Normal 

CD 

Drag 

CL 

Lift 

Cy 

CY 

Side 

X 

cp 

X-C.P. 

Center  of  Pressure  (Cal.  from  MRC) 

Moment  Coefficients: 

CLL 

Roll 

CM 

Pitch 

CLN 

Yaw 

Derivatives: 

^Na 

CNA 

Normal  Force  coef.  deriv.  wrt  Alpha 

mM 

CYB 

Side  Force  coef.  deriv.  wrt  Beta 

^ma 

CMA 

Pitching  Moment  coef.  deriv.  wrt  Alpha 

CLNB 

Yawing  Moment  coef.  deriv.  wrt  Beta 

Cip 

CLLB 

Rolling  Moment  coef.  deriv.  wrt  Beta 

Dvnamic  Derivatives: 

C"/ 

Ip 

Roll  Moment  coef.  deriv  wrt  to  Roll  Rate 

^Nq 

CNQ 

Normal  Force  coef.  deriv  due  to  Pitch  Rate 

^Na 

CNAD 

Normal  Force  coef.  deriv  wrt  to  accel  in  AOA 

C 

mq 

CMQ 

Pitching  Moment  coef.  deriv.  wrt  Pitch  Rate 

c  . 

ma 

CMAD 

Pitching  Moment  coef.  deriv.  wrt  accel.  in  AOA 

2.3  PA  SYSTEM 

The  PA  system  is  also  a  body-axis  system  but  the  orientation  of  the  vehicle  relative  to  the 
resultant  wind,  [Vb^],  is  defined  by  the  total  angle  of  attack,  a'  and  aerodynamic  roll  angle,  (j)',  as 
shown  in  Figure  2.  The  total  angle  of  attack  is  defined  by  the  angle  between  the  longitudinal  axis 
of  missile  and  resultant  wind.  The  total  angle  of  attack  is  always  positive.  The  aerodynamic  roll 
angle  is  defined  as  the  angle  between  the  reference  plane  and  plane  of  resultant  wind  and  is 
considered  positive  if  the  rotation  about  the  body-fixed  x-axis  is  positive  (clockwise  when  looking 
in  the  +X  direction),  thus  a  positive  is  shown. 

The  following  equations  convert  between  the  CA  and  PA  orientation  angles. 

a  =  tan  -l  (tan  a'  cos  <j)') 
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P  =  sin  *1  (sin  a'  sin  (j)') 
a'  =  cos  (cos  a  cos  P) 

(])'  =  tan  (tan  P  /  sin  a) 

The  PA  6  DOF  aerodynamic  equations  for  a  missile  with  tetragonal  symmetry  are  as  follows.'^ 
Force  Coefficients 

C^=Q„(M)  +  C.  (M)a'+C  (M)8^  ;  8^=1^^ 

c;  =  ACy (M,  a')  sin  4(\>'  +  Cy  (M)8r 

‘>9 

C^=C^^(M,a')  +  AC',(M,a')&in^2^'+C'  iM)6q 
Moment  Coefficients 

C,  (M)a'"sin4(|i'+C,^(M)^+C,^^(M)8p 

c  =  c;  (M,a') + AC;,,  (M.aOsin^  2(1.'  +  c;  + C'  (M)8? 

?  2y  °9 

C'  =  Aq,.(W,a')sin4.|.'+C;(M):f^+C:  (M)8r 

The  force  and  moment  coefficients  used  in  the  PA  equations  above  are  listed  in  Table  2  below, 
along  with  the  variable  name  used  in  the  CTA  code. 

2.3  CONTROL  SURFACE  DEFLECTIONS 

The  MD  runs  used  to  calculate  the  coefficients  in  Table  2  require  roll,  pitch  and  yaw  control 
surface  deflections.  Care  must  be  taken  to  insure  that  the  proper  deflections  are  specified  because 
fin  numbering  and  definition  of  a  positive  deflection  varies.  In  MD,  a  positive  deflection  angle 
produces  a  negative  body  axis  rolling  moment  at  zero  angle  of  attack.  The  roll,  pitch,  and  yaw 
control  for  the  MD  convention  can  be  defined  in  terms  of  individual  surface  deflections  by  the 
formulas  for  5p,  5q,  8r  given  below  and  illustrated  in  Figure  3. 

g  —  ~  ^2  ~  ^3  ~  ^4 

P  4 

g  _  -5j  -62+53+64 
9  4 

g  —  +  ^2  +  63  —  64 

^  4 
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Table  2.  PA  Coefficients  and  CTA  Variable  Names. 


Axial  Force  Coefficient: 

Definition: 

CA 

Axial  Force  coef.  at  zero  a’  -  f{M#} 

a' 

CAA 

Variation  of  with  a'  -  f{M#} 

Ca 

h'^eff 

CAD 

Variation  of  with  Effective  Control  Surface  Deflection  -  f{M#} 

Side  Force  Coefficient 

CYP 

Side  Force  coef.  at  =  22.5°  -  f{M#,  a'} 

Cy 

hr 

=  CNr)Q 

Variation  of  Cy  with  Yaw  Control  -  f{M#} 

Normal  Force  Coefficients: 

CN 

Normal  Force  coef.  at  a'  -  f{M#,  a'} 

CMP 

Variation  of  C;  with  roll  angle  (C^at  <t)'=45°-C;at  (t)'=0°)  -  f{M#,  a'} 

c 

5q 

CNDQ 

Pitch  Control  Effectiveness  -  f{M#} 

Rollins  Moment  Coefficients: 

CLLAP 

Induced  Roll  Moment  -  f{M#} 

Cl 

tp 

CLLP 

Roll  Damping  -  f[M#} 

Cl 

%p 

CLLDP 

Roll  Control  Effectiveness  ’  f{M#} 

Pitching  Moment  Coefficients: 

c 

^  m 

CLM 

Pitching  Moment  coef.,  wind  axis  -  f{M#,  a'} 

W,(p 

CLMP 

Variaion  of  with  roll  angle  at  (1)’=45°  -  at  <t)'=0')  ’  f{M#,  a’} 

C'm 

ftiq 

CLMQ 

Pitch  Damping  ‘  f{M#} 

C 

5q 

CLMDQ 

Pitch  Control  Effectiveness '  f{M#} 

Yawing  Moment  Coefficients: 

CLNP 

Yawing  Moment  coef  at  (j)'=22.5°  ’  f{M#,  a'} 

c 

=  CLMQ 

Yaw  Damping '  f{M#} 

c 

"8r 

=  CLMDQ 

Yaw  Control  Effectiveness  ‘  f{M#,} 

3.  CTA  PROGRAM 

The  CTA  program  calculates  the  PA  coefficients  listed  in  Table  2  from  a  standard  MD  output 
file,  FOR006.DAT.  The  program  carries  out  following  steps: 

-  reads  in  data  from  MD  output  file  (FOR006.DAT)  and  the  CTA.DAT  file 

-  outputs  MD  data  to  file  (CART.  OUT)  to  check  that  the  data  has  been  properly  read 

-  calculates  PA  coefficients  from  C A  coefficients 

-  writes  PA  coefficients  to  output  file  (AEROB.OUT) 

A  detailed  explanation  and  instructions  are  in  the  following  sections. 
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3.1  MD 

The  MD  semi-empirical  aerodynamic  prediction  code  estimates  the  aerodynamics  of  a  wide 
variety  of  missile  configurations  to  an  accuracy  suitable  for  preliminary  missile  design^.  Revision 
6/93  is  used  in  this  project. 

Before  running  MD,  an  input  file  (FOR005.DAT)  must  be  created.  It  includes  geometry  and 
flight  conditions.  The  FOROOS.DAT  file  used  in  this  study  is  located  in  Appendix  A  with  a  partial 
listing  of  the  resulting  output  file  (FOR006.DAT)  in  Appendix  B.  The  general  information 
required  to  run  MD  is  covered  in  Reference  5,  so  only  the  specific  requirements  for  CTA  will  be 
covered  here. 

The  angle  of  attack  sweep  and  the  Mach  number  for  each  dataset  are  in  the  FLTCON  section 
of  FOROOS.DAT.  Looking  at  Appendix  A,  note  that  both  the  NMACH  and  NALPHA  inputs  are 
two  digits  (8  is  input  as  08).  This  two  digit  format  must  be  retained  for  these  two  values  to  be 
properly  read  by  CTA  (no  other  variables  have  this  requirement).  There  are  two  ways  to  define 
the  orientation  angles  within  MD;  ALPHA  and  BETA  or  ALPHA  and  PHI.  If  PHI  is  input  and 
non- zero,  it  is  assumed  that  ALPHA  is  the  total  angle  of  attack  (a')  and  PHI  is  the  aerodynamic 
roll  angle  (([)').  The  CTA  code  requires  the  second  method. 

The  baseline  case  is  a  roll  angle  ((})')  of  zero  (fins  are  located  45°  from  the  Y-Z  axis  at  (j)'  =  0). 
The  next  two  cases  set  <[)'  to  22.5°  and  45°  respectively.  Next,  (j)'  is  returned  to  zero  and  the  four 
fin  deflection  cases  are  run,  with  5p,  5q,  5i-  and  6eff  set  to  5°  respectively.  For  the  last  case,  5p  is 
set  to  5°  and  ([)'  to  22.5°.  This  order  must  be  retained.  The  total  number  of  data  sets  (NDS),  64 
in  this  case,  will  change  only  if  the  number  of  Mach  numbers  change  (NDS  =  8  *  NMACH). 
Table  3  below  illustrates  the  required  sequence. 


Tables.  MD  Run  Sequence. 


DATASET 

<t>’ 

5 

1  8  (8  M#) 

0 

0 

9^16 

22.5 

0 

17^24 

45 

0 

25^32 

0 

5p  =  5° 

33-^40 

0 

5q  =  5° 

41  ^48 

0 

5r  =  5° 

49^56 

0 

II 

o 

57-^64 

22.5 

5p  =  5° 

3.2  CTA 

The  next  step  is  to  edit  the  CTA.DAT  file  (Appendix  C).  It  contains  the  total  number  of 
datasets  (NDS),  the  distance  from  the  centerline  to  centroid  of  control  fin  (Ycent)’  the 
number  of  the  fin  set  that  is  being  deflected  (NFS).  For  this  case,  NDS  is  64,  Y^ent  is  0.3751  ft 
(must  be  same  units  as  l^gf)  and  NFS  is  2. 

Finally,  the  CTA.FOR  program  is  compiled,  linked  and  executed.  The  two  input  files  required 
when  rurming  CTA.FOR  are  CTA.DAT  and  FOR006.DAT  previously  discussed.  When  CTA  has 
executed  properly,  two  output  files  will  be  created,  CART. OUT  and  AEROB.OUT.  CART.OUT 
(Appendix  D)  contains  the  data  read  in  from  the  MD  output  file  in  a  headerless  format  to  verify 
input  from  FOR006.DAT.  AEROB.OUT  (Appendix  E)  contains  the  PA  coefficients  in  a  format 
similar  to  that  required  by  the  CAD  AC  flyout  simulation.  The  changes  that  have  to  be  made  are 
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adding  a  slash  at  the  end  of  each  dataset  and  shortening  some  lines  to  less  than  72  columns.  The 
definitions  of  the  variable  names  used  are  listed  in  Table  2. 

All  data  generated  by  MD  and  CTA  should  be  reviewed  for  reasonableness.  If  the 
AEROB.DAT  data  looks  questionable,  the  first  thing  that  should  be  checked  is  a  match  between 
FOR006.DAT  and  CART.OUT.  If  differences  are  found,  CTA  has  not  read  the  input  data 
properly.  This  program  is  not  robust  so  it  is  important  that  the  user  verify  that  the  data  is  being 
read  in  properly.  This  is  especially  important  if  versions  of  MD  other  than  6/93  are  used. 

3.3  CTA  METHODOLOGY 

The  CTA  program  reads  data  from  the  MD  output  file  (FOR006.DAT)  and  inserts  it  into  two 
arrays.  A  one  dimensional  array,  ALPHA,  contains  the  angle  of  attack  sweep  data  and  the  three 
dimensional  BC  array  contains  all  other  input  data.  The  first  and  second  dimension  of  BC  (row 
and  column)  are  listed  in  the  table  below.  The  first  column  contains  general  information  about 
that  dataset  while  the  next  15  columns  contain  the  coefficients  over  the  AOA  sweep  contained  in 
the  ALPHA  array.  The  third  dimension  of  BC  is  the  dataset  number.  Up  to  20  angles  of  attack, 
16  Mach  numbers  and  100  data  sets  can  be  handled  with  the  current  program  dimensions. 

Table  4  BC  Array. 


MD  calculates  all  the  coefficients  listed  in  Table  1  except  so  an  alternate  method  of 

calculating  this  coefficient  is  needed.  A  method  which  predicts  the  roll  damping  derivatives  of 
cruciform-tailed  missiles  was  foimd  in  reference  6,  using  the  following  equation; 

Cip  =  -2.15(rc£:iVT/ Iref)Ci^^ 
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It  is  based  on  empirical  correlation  of  experimental  data  for  several  cruciform-tailed  missiles  at 
Mach  numbers  from  0  to  4.0.  Y^ent’  the  radial  distance  from  body  centerline  to  centroid  of  area 
of  exposed  tail  fin,  is  input  in  the  CTA.DAT  file. 


4.  DATA  COMPARISONS 

This  section  will  compare  the  CTA  generated  PA  aerodynamic  coefficients  to  those  found  in 
the  Rockwell  reports  (Ref.  2  &  Ref.  3).  Rockwell  used  the  ALSAC  computer  code  and  Rockwell 
International  Missile  Drag  and  Aerodynamic  Manuals  to  generate  their  aerodynamic  data.  Wind 
tunnel  data  in  both  PA  and  C A  format  would  be  the  best  way  to  verify  that  the  conversion  process 
in  the  CTA  code  is  working  properly  since  differences  can  be  attributed  to  the  various  methods 
used  by  Rockwell  and  by  MD. 

Figure  4  compares  versus  angle  of  attack  at  two  Mach  numbers.  The  CTA  data  has  a 
smaller  slope  than  the  Rockwell  data  and  shows  little  variation  with  Mach  number  at  angles  of 
attack  less  than  14°. 

Figure  5  compares  versus  angle  of  attack  at  two  Mach  numbers.  Rockwell  predicts 

higher  values  over  the  angle  of  attack  range  at  the  lower  Mach  number  but  the  data  comes  much 
closer  to  matching  at  the  higher  Mach  number. 

Figure  6  compares  versus  angle  of  attack  at  two  Mach  numbers.  The  Rockwell  data 
shows  very  little  variation  with  Mach  number,  the  CTA  data  only  slightly  more.  The  slope  trends 
are  also  different  but  the  overall  agreement  is  quite  good. 

Figure  7  shows  versus  angle  of  attack  at  two  Mach  numbers.  The  slope  of  the  CTA 

data  decreases  at  a  much  higher  rate  with  increasing  angle  of  attack  than  does  the  slope  of  the 
Rockwell  data. 

Figure  8  shows  AC^^,  versus  angle  of  attack.  Although  the  overall  trends  are  similar,  the 

Rockwell  method  predicts  larger  magnitude  values  at  the  low  Mach  number  and  smaller 
magnitude  values  at  high  Mach  number  than  does  the  CTA  method. 

Figure  9  compares  AC'  versus  angle  of  attack  at  two  Mach  numbers.  CTA  predicts  higher 

values  as  angle  of  attack  increases  and  shows  a  larger  variation  with  Mach  number  than  the 
Rockwell  data. 

Figure  10  shows  versus  Mach  number.  The  data  matches  closely  in  the  transonic  region 
but  diverges  in  the  subsonic  and  supersonic  region. 

Figure  11  shows  ^  versus  Mach  number.  Although  the  data  fluctuates  about  the  same 

region,  the  trends  are  similar  only  at  the  higher  Mach  numbers. 

Figure  12  compares  C^  versus  Mach  number.  The  Rockwell  data  shows  no  variation  with 

Mach  number  while  CTA  predicts  a  curve  similar  in  shape  to  the  C^^  versus  Mach  number  curve 
found  in  Figure  10,  with  a  transonic  peak. 

Figure  13  shows  versus  Mach  number.  The  extremes  of  Mach  number  have  similar 

Sq 

values  but  the  transonic  region  is  predicted  to  be  greater  with  the  CTA  method. 

Figure  14  shows  C'  versus  Mach  number.  The  results  are  similar  to  those  in  Figure  12,  with 

dq 

the  CTA  transonic  peak  of  greater  magnitude. 
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Figure  15  shows  versus  Mach  number.  The  supersonic  trends  and  magnitudes  match 

very  well  but  the  transonic  magnitudes  differ  by  a  substantial  margin. 

Figure  16  shows  versus  Mach  number.  The  overall  trends  are  similar  but  the 

magnitudes  are  not. 

Figure  17  compares  C,  versus  Mach  number.  Again,  the  trends  are  similar  but  the 

bp 

magnitudes  differ  substantially. 

Figure  18  shows  versus  Mach  number.  The  magnitudes  differ  substantially  but  are  similar 
in  that  both  changing  little  across  the  Mach  number  range. 


5.  CONCLUSIONS 

The  data  comparison  shows  definite  differences  in  results  although  general  trends  were  often 
similar.  These  differences  were  expected  given  that  both  data  are  estimations  utilizing  different 
methods.  For  a  true  test  of  CTA,  "verified"  aerodynamic  data  in  both  CA  and  PA  format  is 
required  since  CTA  only  converts  MD  data  to  a  different  format. 

The  ability  to  quickly  and  accurately  convert  MD  aerodynamic  data  to  the  PA  format  is 
important  if  PA  equations  are  used  in  a  flyout  simulation.  Although  more  testing  is  in  order,  the 
CTA  code  appears  to  accomplish  the  conversion  process. 
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Figure  4.  Q  versus  Angle  of  Attack. 
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Figure  5.  AC^  vereus  Angle  of  Attack. 
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Figure  6.  versus  Angle  of  Attack. 
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Figure  7.  AC^  versus  Angle  of  Attack. 
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Figure  11.  C.  versus  Mach  Number, 
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Figure  12.  vereus  Mach  Number. 


Figure  13.  Cl  versus  Mach  Number. 
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Figure  16.  C^  .,  versus  Mach  Number. 
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Figure  17.  C,  versus  Mach  Number. 
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APPENDIX  A  -  FOR005.DAT 


CASEID  ROCKPRN 

$FLTCON  NMACH-08.,  MACH-0.2,0.6,0.85,0.95,1.0,1.1,1.6,2.0, 
NALPHA-11.,  ALPHA=0.,2.,4.,6.,8.,10.,12.,14.,16.,18 
ALT-0000 . 0 , PHI-0 . 0 , $ 

$REFQ  XCG-5.4583,  ZCG-0.,  LREF-1.25,  SREF-1.23, 
BLAYER-NATURAL,  RHR-125.,  SCALE-1. 0,$ 

$AXIBOD 


* 


$FINSET1 


$FINSET2 


$DEFLCT 


TNOSE 

= 

CONICAL, 

LNOSE 

- 

2.0917, 

DNOSE 

1.0833, 

BNOSE 

- 

0.375, 

XO 

- 

0.0, 

LCENTR 

- 

8.5, 

DCENTR 

- 

1.25, 

TAFT 

- 

CONICAL, 

LAFT 

- 

1.4483, 

DAFT 

- 

0.6021, 

DEXIT 

*= 

0.0,$ 

SECTYP 

at 

HEX, 

ZUPPER 

- 

0.008,0.0183,0.0432, 

LMAXU 

- 

0.7,0.3,0.25, 

LFLATU 

S 

0.0, 0.0, 0.0, 

SSPAN 

- 

0.0,  0.5833,1.1667, 

CHORD 

- 

3.9167,  1.2833,0.3667, 

XLE 

as 

4.0155, 

SWEEP 

= 

78.5,  57.4, 

STA 

0.,  o.,o.. 

NPANEL 

= 

4., 

PHIF 

8B 

45.,  135.,  225.,  315.,$ 

SECTYP 

HEX, 

ZUPPER 

= 

0.0469,0.0156, 

LMAXU 

= 

0.33,0.25, 

LFLATU 

- 

0.0, 0.0, 

SSPAN 

* 

0.,  0.9167, 

CHORD 

CB 

0.6667,  0.3333, 

XLE 

= 

11.1917,11.525, 

SWEEP 

= 

• 

o 

• 

o 

STA 

ae 

0 . ,  1 . , 

NPANEL 

- 

4., 

PHIF 

= 

45.,  135.,  225.,  315.,$ 

DELTA2 

0.,  0.,  0.,  0., 

XHINGE 

- 

6.7209,11.4584,$ 

DIM  FT 
DERIV  DEG 
SOSE 
DAMP 
SAVE 

NEXT  CASE 

CASEID  ROCKII  -  PHI  22.5 
$FLTCON  PHI-22. 5, $ 
DERIV  DEG 
DAMP 


,20. 
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SAVE 

NEXT  CASE 

CASEID  ROCKII  -  PHI  45 
$FLTCON  PHI-45. 0,$ 

DERIV  DEG 

DAMP 

SAVE 

NEXT  CASE 

CASEID  ROCKII  -  dp  =  5  deg 
$FLTCON  PHI-0. 0,$ 

$DEFLCT  DELTA2  =  -5.,  -5.,  -5.,  -5.,$ 

DERIV  DEG 
DAMP 
SAVE 

NEXT  CASE 

CASEID  ROCKII  -  dq  =  5  deg 
$DEFLCT  DELTA2  =  -5.,  -5.,  5.,  5.,$ 

DERIV  DEG 
DAMP 
SAVE 

NEXT  CASE 

CASEID  ROCKII  -  dr  -  5  deg 
$DEFLCT  DELTA2  -  -5.,  5.,  5.,  -5.,$ 

DERIV  DEG 
DAMP 
SAVE 

NEXT  CASE 

CASEID  ROCKII  -  deff  =  5  deg 
$DEFLCT  DELTA2  =-5.,  -5.,  5.,  5.,$ 

DERIV  DEG 
DAMP 
SAVE 

NEXT  CASE 

CASEID  ROCKII  -  dp  =  5  deg,  phi  =22.5  deg 
$FLTCON  PHI-22. 5, $ 

$DEFLCT  DELTA2  -  -5.,  -5.,  -5.,  -5.,$ 

DERIV  DEG 
DAMP 

NEXT  CASE 
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APPENDIX  C  -  CTA.DAT 
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0.3751 

2 


-  #  OF  DATA  SETS  -  NDS 

-  DISTANCE  FROM  CENTERLINE  TO  CENTROID  OF  FIN  (Iref  units)  -  YCENT 

-  NUMBER  OF  FIN  SETS  -  NFS 
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.141E+07  0.00  0.00  1.230  1.250  5.458 
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APPENDIX  F  -  CTA.FOR 


C  CTA4.F0R  29  NOV  94 

C  THIS  PROGRAM  CONVERTS  MISSILE  DATCOM  OUTPUT  IN  CARTESIAN  COORDINATES 
C  TO  AEROBALLISTIC  COORDINATES.  GREG  WILDER 
C  ASC/XREWS  904-882-3722 

C23456789012345678921234567893123456789412345678951234567896123456789712 
C 

DIMENSION  ALPHA(20),BC(16,20,100),CN(16,20),CNP(16,20), 

&  CLM(16,20),CLMP(16,20),CYP(16,20),CLNP(16,20),CA(16), 

&  CAA( 16 ) ,CAD( 16 ) , CNDQ( 16 ) ,CLMDQ( 16 ) ,CLLDP( 16 ) , CLLP( 16 ) , 

&  CLMQ(16),CLLAP(16) 

REAL  LREF 

CHARACTER*6  ANUM,ADYN,ACASE 
C 

C  OPEN  INPUT  FILE  &  OUTPUT  FILE 

OPEN(UNIT=5,FILE='CTA.DAT' ,STATUS='OLD' ) 

OPEN( UNIT-6, FILE- 'FOR006.DAT' ,STATUS='OLD' ) 
OPEN(UNIT-9,FILE='CART.OUT' ,STATUS='NEW' ) 
OPEN(UNIT-10,FILE-'AEROB.OUT' ,STATUS='NEW' ) 

C 

C  VARIABLES: 

C  ALPHA  -  ALPHA  MATRIX 

C  BC  -  BODY  COEFFICEINTS  FROM  DATCOM  -  MAX  20  ALPHA,  MAX  16  M# 

C  &  100  DATA  SETS 
C  NDS  -  #  OF  DATA  SETS 

C  YCENT  -  DISTANCE  FROM  BODY  CENTERLINE  TO  CENTROID  OF  CONTROL  FIN 
C  SAME  UNITS  AS  LREF 

C  NFIN  -  NUMBER  OF  THE  FIN  SET  TO  BE  DEFLECTED 
C  RAD  -  DEGREES  TO  RADIANS  CONVERSION  FACTOR 
C  DEG  -  RAD  TO  DEG  CONVERSION  FACTOR 

C  122  -  FACTOR  REQUIRED  TO  SKIP  TO  RUN  WITH  EQUAL  M#  @22.5  DEG 
C  145  -  FACTOR  REQUIRED  TO  SKIP  TO  RUN  WITH  EQUAL  M#  @45  DEG 
C 

RAD  «  1.0/57.2958 
DEG  =  57.2958 
READ(5,*)NDS 
READ (5,*) YCENT 
READ(5,*)NFIN 
C 

C  READ  IN  #  OF  M#  &  #  OF  AOA'S  (LINE  16  &  17)  FROM  INPUT  DECK  ECHO 
C 

DO  100  ILOOK-1,999 
READ(6,1000)ACASE 
1000  FORMAT(5X,A6) 

100  IF  (ACASE  .EQ.  'CASEID' )GOTO  10 
C 

C  READ  IN  NMACH  &  NALP  (REAL  IN  BC) 

C 

10  READ(6,1010)BC(2,1,1) 

WRITE(9,*)BC(2,1,1) 

1010  FORMAT(21X,F2.0) 

READ( 6, 1015 )BC( 1,1,1) 

WRITE(9,*)BC(1,1,1) 

1015  FORMAT(23X,F2.0) 

C 

C  BEGIN  READ  OF  DATA 
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c 

DO  110  I-1,NDS 
C 

C  #  OF  M#  &  AOA'S  ARE  COJSTANT 

BC(1,1,I)-BC(1,1,1) 

BC(2,1,I)=BC(2,1,1) 

C 

C  LOCATE  FIRST/NEXT  DATA  SET  BY  FINDING  "NUMBER"  (IN  MACH  NUMBER) 

C 

DO  120  ILOOK«=l,  10000 
15  READ(6,1020)ANUM 
1020  FORMAT(lX,A6) 

IF(ANUM  .EQ.  'NUMBER' )C30T0  20 
120  CONTINUE 
C 

C  READ  HEADER  INFO  AND  INSERT  INTO  FIRST  COL  OF  BC  MATRIX 
C 

20  READ(6,1040)BC(3,1,I),BC(4,1,I),BC(5,1,I),BC(6,1,I),BC(7,1,I), 

&  BC(8,1,I),BC(9,1,I),BC(10,1,I),BC(11,1,I) 

C 

WRITE(9,1041)BC(3,1,I)/BC(4,1,I),BC(5,1,I)/BC(6,1/I),BC(7,1,I), 

&  BC(8,1,I),BC(9,1,I),BC(10,1,I),BC(11,1,I) 

C 

1040  FORMAT(/,2X,F4.2,2X,F10.2,37X,E9.3,2X,F8.2,2X,F7.2,2X,F8.3,2X, 

&  F7.3,11X,F8.3,2X,F9.3,//////) 

1041  FORMAT(/,2X,F4.2,2X,F10.2,37X,E9.3,2X,F8.2,2X,F7.2,2X,F8.3,2X, 

&  F7.3,11X,F8.3,2X,F9.3,/) 

C 

C  READ  FIRST  BLOCK  OF  DATA 

C23456789012345678921234567893123456789412345678951234567896123456789712 
NALPHA-IFIX  (  BC( 1,1,1)  ) 

DO  125  J-1,NALPHA 

READ(6,1050)ALPHA(J),BC(J,2,I),BC(J,3,I),BC(J,4,I),BC(J,5,I), 

&  BC(J,6,I),BC(J,7,I),BC(J,8,1),BC(J,9,I),BC(J,10,I),BC(J,11,I), 

&  BC(J,12,I) 

WRITE(9,1050)ALPHA(J),BC(J,2,I),BC(J,3,I),BC(J,4,I), 

&  BC(J,5,I),BC(J,6,I),BC(J,7,I),BC(J,8,I),BC(J,9,I), 

&  BC(J,10,I),BC(J,11,I),BC(J,12,I) 

1050  FORMAT(3X,F5.2,1X,3(2X,F7.3),4X,3(2X,F7.3),4X,5(2X,E10.3)) 

C 

125  CONTINUE 
C  READ  IN  XCP 

READ(6,1055) 

1055  FORMAT(//) 

DO  130  J=1,NALPHA 
READ(6,1060)BC(J,13,I) 

WRITE(9,1061)BC(J,13,I) 

1060  FORMAT(78X,F8.3) 

1061  FORMAT(8X,F8.3) 

C 

130  CONTINUE 

C 

C  READ  IN  FIN  DEFLECTIONS  (IF  NFIN  =  2  SKIP  TO  SECOND  FIN  SET) 

C 

READ(6,1065) 

1065  FORMAT(/) 

C  IF(NFIN  .EQ.  2)READ(6,1065) 

132  READ(6,1066)IFS,BC(12,1,I),BC(13,1,I),BC(14,1,I),BC(15,1,I) 

IF  (IFS  .LT.  NFIN)GOTO  132 

WRITE(9,1067)IFS,BC(12,1,I),BC(13,1,I),BC(14,1,I),BC(15,1,I) 
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1066  F0RMAT(43X,I2,4(4X,F7.2)) 

1067  F0RMAT(3X,I2,4(4X,F7.2)) 

C 

C  READ  IN  DYNAMIC  DERIV 

C  -IXXATE  DYNAMIC  DERIV  DATA  BY  FINDING  "DYNAMI"  (IN  HEADER) 

C 

DO  140  ILOOK=l, 10000 
READ(6,1070)ADYN 
1070  FORMAT(49X,A6) 

IF(ADYN  .EQ.  'DYNAMI' )G0T0  30 
140  CONTINUE 
C 

30  READ(6,1075) 

1075  FORMAT(/) 

DO  150  J=1,NALPHA 

READ(6,1080)BC(J,14,I),BC(J,15,I),BC(J,16,I) 

WRITE(9,1085)BC(J,14,I),BC(J,15,I),BC(J,16,I) 

1080  FORMAT(47X,E11.3,10X,E10.3,11X,E13.5) 

1085  FORMAT(8X,E10.3,10X,E10.3,11X,E13.5) 

150  CONTINUE 
C 

110  CONTINUE 
C 

C  CALCULATE  CN  -  AEROBALLISTIC  CN  PRIME 

C  CNP  -  CORRECTION  TO  CN  WHEN  PHIP  IS  NOT  ZERO.  DELTA  CNP 

C  *  CN  @  PHIP  =  45  -  CN  @  PHIP  =  0 

C  TCNP45  -  TEMP  CN  @  PHIP  =  45  DEG 
C  CLM  -  AEROBALLISTIC  CM  -  @  PHIP  =  0 

C  CLMP  -  CORR  TO  CLM  WHEN  PHIP  IS  NOT  ZERO. 

C  DELTA  CLM  =  CLM  @  PHIP  =  45  -  CLM  @  PHIP  =  0 

C  CA  -  CA  -  AXIAL  COEF  AT  ZERO  AOA 

C 

NMACH=IFIX(BC(2,1,1)) 

WRITE( 9 , 1100 )NMACH,BC( 2 , 1 , 1 ) 

1100  FORMAT(2X, 'NMACH  «  ' ,14, 'BC(2,1,1)  =  ',F10.4) 

122  =  NMACH 
145  -  2*NMACH 
C 

DO  200  1=1, NMACH 
DO  210  J=1,NALPHA 
C 

CN(I,J)  =  BC(J,2,I) 

C 

PHIP  =  RAD  *  BC(7, 1,1+145) 

TCNP45  =  BC(J, 2,1+145)  *  COS(PHIP)  -  BC(J, 5,1+145)  *  SIN(PHIP) 
CNP(I,J)  =  TCNP45  -  BC(J,2,I) 

C 

CLM(I,J)  =  BC(J,3,I) 

TCMP45  =  BC(J, 3,1+145)  *  COS(PHIP)  -  BC(J, 6,1+145)  *  SIN(PHIP) 
CLMP(I,J)  =  TCMP45  -  BC(J,3,I) 

C 

PHIP  =  RAD  *  BC(7, 1,1+122) 

C 

CYP(I,J)  =  BC(J, 5,1+122)  *  COS(PHIP)  + 

&  BC(J, 2,1+122)  *  SIN(PHIP) 

CLNP(I,J)  =  BC(J,3,I+I22)*SIN(PHIP)+BC(J,6,I+I22)*C0S(PHIP) 

C 

210  CONTINUE 
C 

CA(I)  =  BC(1,4,I) 
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c 

200  CCOTINUE 
C 

C  CALC  CAA  &  CAD 
C 

DO  300  I«1,NMACH 
C 

C  CAA 
C 

CASUM  =  0.0 
DO  310  J-1,NALPHA-1 

DALP=  ALPHA(J+1)-ALPHA(J) 

CASUM  -  (BC(J+1,4,I)  -  BC(J,4,I))/DALP  +  CASUM 
310  CCOTINUE 
C 

CAA(I)  =  CASUM  /  (NALPHA-1) 

C 

C  CALC  DQ,  DR  &  DEFF 
C 

D1  -  BC(12,1,I+6*NMACH) 

D2  -  BC(13,1,I+6*NMACH) 

D3  -  BC( 14,1, 1+6 *NMACH) 

D4  -=  BC(15,1,I+6*NMACH) 

DQ  =  (-Dl-D2+D3+D4)/4.0 
DR  =  (-Dl+D2+D3-D4)/4.0 
DEFF  >=  (  ABS(DQ)  +  ABS(DR)  )  /  2.0 
C 

C  CAD  (CALC  @  ZERO  AOA  ONLY. ) 

C 

J  -  1 

CAESUM  -  (BC(J,4,I+6*NMACH)  -  BC(J,4,I)) 

C 

CAD(I)  -  CAESUM  /  (DEFF**2) 

C 

C  CNDQ  &  CLME)Q  LOOP  @  ZERO  AOA 
C 

C  CALC  DQ 
C 

Dl  »  BC( 12,1, 1+4 *NMACH) 

D2  »  BC(13,1,I+4*NMACH) 

D3  -=  BCU4,1,I+4*NMACH) 

D4  -  BC(15,1,I+4*NMACH) 

DQ  *  (-Dl-D2+D3+D4)/4.0 
C 

J-1 

CNDQSUM  »  (BC(J,2,I)  -  BC( J,2,I+4*NMACH) ) 

CLMDQSUM  =  (BC(J,3,I)  -  BC( J, 3,1+4 *NMACH) ) 

C 

CNDQ(I)  =  CNDQSUM/(DQ) 

CLMDQ(I)  -  CLMDQSUV(DQ) 

C 

C  CLLAP 
C 

CLSUM  =  0.0 
NSUM  =  0 

DO  325  J=1,NALPHA 
ACLL  -=  ABS(BC(J,7,I+NMACH)) 

IF  (  ACLL  .LT.  0.001  )  GO  TO  325 
C  ONLY  NCW-ZERO  CLL  VALUES  ARE  INCLUDED  IN  THIS  CALCULATION 

CLSUM  «  ((BC(J,7,I+NMACH)-BC{J,7,I))/ALPHA{J)**2)  +  CLSUM 
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NSUM  =  NSUM  +  1 
325  CaWINUE 

CLLAP(I)  =  CLSUM/NSUM 
C 

C  CLLDP  @  ZERO  AQA 
C 

C  CALC  DP 
C 

Dl  =  BC(12,1,I+3*NMACH) 

D2  =  BC(13,1,I+3*NMACH) 

D3  =  BC(14,1,I+3*NMACH) 

D4  =  BC(15,1/I+3*NMACH) 

DP  -=  (-Dl-D2-D3-D4)/4.0 
C 

J>=1 

CLLDP(I)  =  (BC(J,7,I+3*NMACH)  -  BC( J,7, I ) )/(DP) 
C 

C  CALC  CLLP 
C 

LREF  =  BC(9,1,1) 

CLLP(I)  =  -2.15*(YCENT/LREF)*CLLDP(I) 

C 

C  CALC  CLMQ 
C 

J=1 

CLMQ(I)  -  BC(J,16,I) 

C 

300  CONTINUE 
C 

C  PRINT  COEFFICIENTS 
C 

C  WRITE  ALPHA  &  M#  INFORMATION 
C 

WRITE( 10,1200 )NALPHA,NMACH 
1200  FORMAT(lX,'C  NALPHA=' ,I3,2X, 'NMACH=' ,I3) 

WRITE(10,1210) (ALPHA( J) , J=1,NALPHA) 

1210  FORMAT(lX,'C  AQA  ACROSS  (DEG) : ' ,20(F4.1, ' , ' ) ) 
WRITE(10,1212) (ALPHA( J) , J=1,NALPHA) 

1212  FORMAT(6X, 'DATA  ALPHA/' f20(F4.1, ) 

WRITE ( 10 , 1220 ) ( BC ( 3 , 1 , J ) , J=1 ,NMACH ) 

1220  FORMAT(lX,'C  M#  DCWN; ' ,16(F4.2, ' , ' ) ) 
WRITE(10,1222)(BC(3,1,J),J=1,NMACH) 

1222  FORMAT(6X,'DATA  MACH/' ,16(F4.2, ' , ' ) ) 

C 

C  WRITE  CN 
C 

WRITE(10,1230) 

1230  FORMAT(7X, 'DATA  TCN/') 

C 

DO  320  I-1,NMACH 

320  WRITE(10,1250)(CN(I,J),J=1,NALPHA) 

C 

C  WRITE  CNP 
C 

WRITE(10,1240) 

1240  FORMAT(  7X, ' DATA  TCNP/' ) 

DO  340  I=1,NMACH 

340  WRITE(10,1250)(CNP(I,J),J=1,NALPHA) 

C 

C  WRITE  CLM 
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c 

WRITE(10,1245) 

1245  F0RMAT(7X, 'DATA  TCLM/' ) 

DO  345  I-1,NMACH 

345  WRITE( 10 , 1250 ) ( CLM( I , J ) , J-1 ,NALPHA) 

C 

C  WRITE  CLMP 
C 

WRITE(10,1260) 

1260  FORMAT( 7X, ' DATA  TCLMP/' ) 

DO  350  I=1,NMACH 

350  WRITE(10,1250)(CLMP(I,J),J=1,NALPHA) 

C 

C  WRITE  CYP 
C 

WRITE(10,1270) 

1270  FORMAT ( 7X , ' DATA  TCYP/' ) 

DO  360  I-1,NMACH 

360  WRITE(10,1250)(C:YP(I,J),J=1,NALPHA) 

C 

C  WRITE  CLMP 
C 

WRITE(10,1280) 

1280  FORMAT( 7X, 'DATA  TCLNP/' ) 

DO  370  I=1,NMACH 

370  WRITE( 10 , 1250 ) ( CLNP( I , J ) , J=1 ,NALPHA) 

C 

C  WRITE  CA 
C 

WRITE( 10,1290) 

1290  FORMAT(7X, 'DATA  TCA/') 

C 

WRITE( 10 , 1250 ) ( CA( I ) , 1=1 ,NMACH ) 

C 

C  WRITE  CAA 
C 

WRITE( 10,1300) 

1300  FORMAT(  7X, ' DATA  TCAA/' ) 

WRITE( 10 , 1310 ) ( CAA( I ) , 1=1 ,NMACH) 

1310  FORMAT(5X, '&',20(F8.5,',')) 

1312  FORMAT(5X,'&',20(F9.4,',')) 

C 

1250  FORMAT(5X,'&',20(F7.3,',')) 

C 

C  WRITE  CAD 
C 

WRITE(10,1320) 

1320  FORMAT( 7X, 'DATA  TCAD/' ) 

WRITE( 10 , 1310 ) ( CAD( I ) , 1=1 ,NMACH) 

C 

C  WRITE  CNDQ 
C 

WRITE(10,1330) 

1330  FORMAT(7X, 'DATA  TCNDQ/') 

WRITE( 10 , 1310 ) { CM)Q( I ) , 1=1 ,NMACH) 

C 

C  WRITE  CLMDQ 
C 

WRITE(10,1340) 

1340  FORMAT(7X,  'DATA  TCLMDQ/' ) 
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WRITE( 10 , 1310 ) ( CLMDQ( I ) , 1=1 ,NMACH) 
C 

C  WRITE  CLMQ 
C 

WRITE(10,1370) 

1370  FORMAT(7X, 'DATA  TCLMQ/') 

WRITE( 10 , 1312 ) ( CLMQ( I ) , 1=1 ,NMACH) 

C 

C 

C  WRITE  CLLAP 
C 

WRITE(10,1345) 

1345  FORMAT(7X, 'DATA  TCLLAP/') 

WRITE{ 10 , 1310 ) ( CLLAP( I ) , 1=1 ,NMACH) 
C 

C  WRITE  CLLDP 
C 

WRITE(10,1350) 

1350  FORMAT(7X, 'DATA  TCLLDP/') 

WRITE( 10 , 1310 ) ( CLLDP( I ) , 1=1 ,NMACH) 
C 

C  WRITE  CLLP 
C 

WRITE(10,1360) 

1360  FORMAT(7X, 'DATA  TCLLP/') 

WRITE( 10 , 1310 ) ( CLLP( I ) , 1=1 ,NMACH) 

C 

STOP 

END 
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APPENDIX  G  -  CTA  EQUATION  SUMMARY 


CA  —  Cj^g  —  ) 


CAA-L^,- ^ 


#a-l  Kr 
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^  Aa 


CAD  =  C, 


dCA  AQ  ^  Q(S,ff^5°.a'=0°)  Q(S,jy=0°,a'^0°) 


^82^  35^^  A5^ 


cyp  =  AC;:,^-  =  Cy(,,22.5»)Cos(<^')+<^%=22.5»)Sin(<l)')  =  /(a',M*) 

CA^  =  C;„=Q  =  /(a',M*) 

CATP  =  AQ^,  =  =  m\  M* ) 

^f=45«)  ^  ^^^'=45“)  “  ^V=45») 

*  ^/  %  jf^  \ 

C^®2  =  c;,,=^-^  =  /(M) 


#a-l 

I 

CLLAP  =  Q^-  =  — 

a  4 


Q«.'=22.S°)  ^/(<»'=0°) 

Aa^ 


(#a-l) 


=nM*) 


CLLDP  =  C,  = 

8p 


Q  =  C/8'  ,oc'^sin4())’ 

CLIP  =  =  -lASiYcEifr  !  ~  /(^  ) 

r  — C 

W=5°,«'="0°)  ^(»=0°,«'=0°) 

bp 

CLM  =  C„  =  Q=m',M*) 

CLMP  =  AC^y  =  C;(^,=450)  =  /(a',M*) 

C(4.'=45°)  =  Q(1.'=45”)  COS0'-C„(^.^45O)  Sin^' 
^m(<ti'=0'’)  “  Q(«t>'=0‘’) 

CLMQ  =  +  Q6.(„=o-)  >  =  ) 

CLNP  =  AC;^,  =  C;(^,=22.5»)  =  /(a'’  ) 

^n(it)'=22.5°)  =  ^m((t>'=22.5°)  <t>^  +  Q((t)'=22,5‘’) 


=  /(M") 
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